gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十一章/eg11_3.m
k=0:0.01:100; for index=1:10000 p=[2 15 27 k(index)+12 k(index)+1]; r=roots(p); if max(real(r))>0 break; end end sprintf('系统临界稳定时K值为:K=%7.4f\n',k(index))